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Abstract 

The clonal dissemination of VanB-type vancomycin-resistant Enterococcus faecium (VREfm) strains in three Swedish hospitals 
between 2007 and 201 1 prompted further analysis to reveal the possible origin and molecular characteristics of the 
outbreak strain. A representative subset of VREfm isolates (n = 18) and vancomycin-susceptible E faecium (VSEfm, n = 2) 
reflecting the spread in time and location was approached by an array of methods including: selective whole genome 
sequencing (WGS; n = 3), multi locus sequence typing (MLST), antimicrobial susceptibility testing, virulence gene profiling, 
identification of mobile genetic elements conferring glycopeptide resistance and their ability to support glycopeptide 
resistance transfer. In addition, a single VREfm strain with an unrelated PFGE pattern collected prior to the outbreak was 
examined by WGS. MLST revealed a predominance of ST192, belonging to a hospital adapted high-risk lineage harbouring 
several known virulence determinants (n>10). The VREfm outbreak strain was resistant to ampicillin, gentamicin, 
ciprofloxacin and vancomycin, and susceptible to teicoplanin. Consistently, a v'an82-subtype as part of Tn 7549/Tn5382 with 
a unique genetic signature was identified in the VREfm outbreak strains. Moreover, Southern blot hybridisation analyses of 
PFGE separated SI nuclease-restricted total DNAs and filter mating experiments showed that vanB2-Jr\l 549/Tr\5382 was 
located in a 70-kb sized rep^/pRUM plasmid readily transferable between E faecium. This plasmid contained an axe-txe toxin- 
antitoxin module associated with stable maintenance. The two clonally related VSEfm harboured a 40 kb repu/pR\jM plasmid 
absent of the 30 kb vanB2-Jr\l 549/Jr\5382 gene complex. Otherwise, these two isolates were similar to the VREfm outbreak 
strain in virulence- and resistance profile. In conclusion, our observations support that the origin of the multicentre outbreak 
was caused by an introduction of vanB2-Jr\l 549/Tr\5382 into a rep^/pRUM plasmid harboured in a pre-existing high-risk E 
faecium ST192 clone. The subsequent dissemination of VREfm to other centres was primarily caused by clonal spread rather 
than plasmid transfer to pre-existing high-risk clones. 
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Introduction 

Enterococci, and Enterococcus faecium in particular, have 
undergone a transition from harmless gut commensals to be a 
leading cause of multidrug resistant hospital infections [1]. E. 
faecium is associated with urinary tract infections, endocarditis, 
infections in indwelling catheters and septicaemia in hospitalised 
patients [2,3]. Notably, a pronounced increase in bacteraemias 
caused by E. faecium in Europe has been reported [4]. The ability 
of hospital adapted lineages of E. faecium to compile antibiotic 
resistance and virulence factors by horizontal gene transfer might 
be attributable to this observation [1]. 

The use of Multi Locus Sequence Typing (MLST) has been 
considered a standard method for global epidemiological surveil- 



lance [5], while pulsed-field gel electrophoresis (PFGE) is the 
preferred method for examination of E. faecium seen in local 
outbreaks. The clonal complex (CC) 17 has shown to pool hospital 
associated E. faecium strains characterised by a high rate of 
recombination, multidrug resistance as well as numerous virulence 
determinants [6-8] . Newer insights pertained with another type of 
population structure analysis, Bayesian Analysis of Population 
Structure (BAPS) of CC17 strains, show divergent origins of 
sequence type (ST) lineages within GC17. GC17 strains could 
largely be divided in two BAPS groups, 2-1 (lineage 78), and 3-3 
(lineages 17 and 18), with the corresponding MLST ancestry 
nodes in parenthesis [9] . 
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A worldwide increased prevalence of acquired resistance to 
commonly used antibiotics is observed in clinical isolates of E. 
faecium. A total of eight gene clusters, vanA,B,D,E,G,L,M,N, 
have been associated with acquired vancomycin resistance in 
enterococci (VRE) [10—13]. VanA VRE is most prevalent globally, 
but VanB-type VRE are predominant in Australia and on the rise 
in many European countries [14-19]. The vanB gene cluster 
include three subtypes, vanBl-3, conferring inducible low- to 
high-level resistance to vancomycin and susceptibility to teicopla- 
nin [20]. The predominant subtype vanB2 is an integral part of an 
Integrative Conjugative Element transposon family Tnl549/ 
Tn3382 supporting transfer of vanB2 [21]. Several plasmid 
replicon types encoding glycopeptide resistance as well as 
stabilising toxin-antitoxin systems have been linked to CC17 
strains [22]. Moreover, a number of putative virulence genes have 
been associated with E. faecium and it is hypothesised that their 
phenotypes might work in concert to promote host colonisation 
and subsequent invasion [2]. 

The prevalence of vancomycin resistant E. faecium (VREfm) in 
Sweden remained low until 2007 when a large hospital associated 
outbreak occurred [15]. Three hospitals in separate counties were 
involved. The outbreak was not declared over until 2011. 

The aim of the present study was to explore the origin of the 
outbreak strain by performing molecular characterisation of 
representative isolates of the outbreak strain and compare them 
with consecutive invasive E. faecium isolates from the same time 
period and location. 

Materials and Methods 

Bacterial isolates 

All cases of vancomycin-resistant enterococci (VRE) are 
mandatorily reported to the Public Health Agency of Sweden 
and collected for resistance and epidemiological typing. In 2007, 
an increasing number of notified VRE-cases were seen in 
Stockholm County, related to a clonal VanB-type E. faecium 
strain. The strain was subsequently reported in two other 
geographically distant counties (Vastmanland and Halland) [15]. 
During the autumn 2008, clonally related isolates of VanB-type 
Efm (n = 17) from separate infection/ colonisation events in these 
three counties were selected for molecular studies. The selected 
isolates occurred early and late during the outbreak period and 
represented isolates with divergent resistance profiles and related 
PFGE subtypes. Four blood isolates from Collection B (see below) 
were also included in the molecular analyses, yielding a total of 2 1 
isolates in Collection A. 

In an attempt to reveal the origin of the E. faecium vanB 
outbreak strain, all consecutive E. faecium blood culture isolates 
from of January 2006 to 31'' of August 2009 (n=191; 
Collection B) diagnosed at the Karolinska University Hospital 
Huddinge where the outbreak was first identified, were collected 
and analysed by PFGE; 2006 (n = 45), 2007 (n = 32), 2008 (n = 7 1) 
and 2009 (n = 43). Four of these isolates were selected for further 
analysis using whole genome sequencing (WGS) with the Roche 
454 pyrosequencing platform. The selection criterion was based 
on PFGE patterns, where the first vancomycin susceptible 
(VSE1036), the first vancomycin resistant (VRE 1044) and the 
most recent vancomycin resistant (VRE 1261) isolate with indis- 
tinguishable or closely related PFGE-patterns to the outbreak 
strain, were chosen. For comparison, one vancomycin resistant 
isolate (VRE0576) from 2006 was chosen because of its divergent 
PFGE pattern. E. faecium 64/3 [23], BM4105RF and BM4105- 
Str [24] were used as recipients in filter mating experiments. 
Isolates from a polyclonal cluster of vanB2 positive E. faecium 



from 2002-2004 in the Swedish county Orebro [18] were included 
in the ICE5/itvan Q8 PGR to evaluate their vanB transposon 
signature. 

Antimicrobial susceptibility testing 

The minimum inhibitory concentration (MIC) was determined 
using Etest (BioMerieux) and interpreted according to the clinical 
breakpoints of the European committee on Antimicrobial Suscep- 
tibility Testing (EUCAST) (www.eucast.org). For ciprofloxacin, a 
tentative breakpoint was used classifying isolates with MIC > 
32 mg/L as high level resistant [25]. 

PFGE and MLST 

For 5m<2l-digestion, the protocol adapted by Saeedi et al. [26] 
was used with 5 U/mL lysozyme added in the lysis buffer. The 
bands were separated with the following program: Block I switch 
time 3 to 26,5s for 14 hours and 50 minutes. Block II: switch time 
0,5 to 8,5s for 6 hours and 25 minutes. Total run time 21 hours 
and 15 minutes at 6V with 120°. The PFGE patterns were 
analysed and compared using BioNumerics software (version 6.6, 
Applied Maths). The Dice coefiicient was used for pair-wise 
comparison of patterns, and the un-weighted pair group method 
with arithmetic mean (UPGMA) for pattern grouping. Isolates 
clustering above 97% were considered identical and isolates with 
identity >90% closely related. 

MLST was performed using the method adapted by Homan et 
al. [5] with the following primers: adkln, adk2n, atpln, atp2n, 
ddll, ddl2, gdhl, gdh2, gydl, gyd2, pstSln, pstS2n, purKln and 
purK2n. 

Detection of genes by PCR and isolation of bacterial DNA 

Extraction of DNA for all PCRs were performed by BioRobot 
M48 (Qiagen), according to the manufacturers manual. Primers 
and positive controls are described in Table S 1 . 

The chosen virulence genes are associated with high-risk 
genotypes, and included esp [27], hyl [28], acm [29], efaAfm 
[30], sgrA, ecbA, scm, orf903/20 10/23 14 and pilA/B [31-33]. 
PCRs genotyping presence of vanB [20] and linkage to Tn3382 
[34] was done. PCRs were conducted as stated in Table SI. PCR 
was performed using the JumpStart REDTaq Readymix PCR 
Reaction mix (Sigma), with a standard program of 1 min in 95°C 
followed by 30 cycles at 95°C for 30 sec, 30 sec of annealing in the 
temperature given in Table SI, 72 °G for 1 min with a final 
elongation step at 72°C for 7 min. The presence of vanB in 
transconjugants was tested by PCR using 1 |il of bacterial culture 
in BHI broth and an additional initial denaturation step of 1 0 min 
at 95°C. 

Southern blotting and hybridisation 

PFGE analysis of SI -digested DNA was used to analyse the 
plasmid content. Plugs were made as for Smal digestion, and the 
digestion was performed as described by RosvoU et al. [22]. The 
Vacugene XL system (Amersham Biosciences) was used for 
Southern blotting. Consecutive hybridisation was performed using 
^^i^i7/pRUM5 vanB and axe-txe probes in the mentioned order. 
Characterisation of the reppi^Qi [35] and Tep2/pKE25 [36] plasmid 
determinants were also done by Southern blotting and hybridisa- 
tion after Sl-nucelase PFGE. Probes were made by amplification 
using positive controls (see Table SI), and labelled using the PCR 
DIG synthesis kit (Boehringer Mannheim). The same hybridisa- 
tion protocol as in RosvoU et al. [22] was used with the following 
modification: The DNA was purified after the first PCR using the 
Cycle Pure Kit (zDNA). 
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Conjugative transfer of vanB 

Filter mating was performed according to Bjorkeng et al. [18] 
with some minor modifications, using the E. faecium 64/3 and E. 
faecium BM4105-RF as recipient strains for the first filter mating, 
and BM4015-Str in retransfer. Briefly, the isolates were grown 
together on MF-Millipore membrane filters for 24 h, spotted on 
selective BHI agar plates containing either vancomycin (8 mg/L), 
fusidic acid (10 mg/L) and rifampicin (20 mg/L), or all three 
antibiotics together. The bacterial suspension was serially diluted 
down to 10~^ and incubated at 37°C for 48 h. In the retransfer 
experiments, the recipients were selected on plates containing 
1000 mg/L streptomycin. 

Whole genome sequencing and analysis 

Chromosomal DNA from the four isolates (VSE1036, 
VRE1044, VRE1261 and VRE0576) was prepared using the 
DNeasy Blood and Tissue Kit (Qjagen) with lysozyme (20 mg/ 
mL) added to the lysis buffer and further treated with RNase. The 
protocol also allowed purification of plasmids. Libraries were 
prepared and used for whole genome shotgun sequencing on a 4- 
region picotiter plate with the Roche 454 FLX system according to 
standard protocols (www.454.com). Raw sequencing data were 
processed with standard filters using the GS Run Processor (v 2.6), 
generating between 246084 and 310421 reads for each of the 
strains with average lengths between approximately 307 and 320 
nucleotides, corresponding to between 78750175 and 99124035 
nucleotides. Reads were assembled de novo with the accompany- 
ing GS de novo assembler software (v 2.6) (454 Newbler algorithm) 
generating between 201 and 302 contigs with a length of more 
than 100 nucleotides. The GS Reference mapper software (v 2.6) 
was subsequently used for homology comparisons between the 
different strains and to find homologies to specific query gene 
sequences and also used to identify indels and point mutations that 
separated the different strains. Some small plasmids could be 
identified by screening for contigs where individual reads mapped 
to both ends of the contig. The contigs were tentatively linked to 
each other by comparison of molecular biology data, identification 
of individual reads mapping to two different contigs and 
comparisons with published genomes. The tentative gene content 
of the contigs was automatically analysed at NGBI with the 
Prokaryotic Genome Annotation Pipeline (PGAP). 

In Figure SI alignment of genome sequences against the 
reference genome of Aus0004 (NC_0 17022) was done by 
reordering the contigs of VRE1044 against the Aus0004 reference 
genome using Mauve 2.3.1. The reordered contigs of VRE1044 
were then used as reference to reorder the contigs of the other 
isolates before the genomes were aligned using progressiveMauve 
[37]. 

The Whole Genome Shotgun projects have been deposited at 
DDBJ/EMBL/GenBank under the accessions JAAJOOOOOOOO 
(VSE1036), JAAKOOOOOOOO (VRE0576), JAALOOOOOOOO 
(VRE1044) and JAAMOOOOOOOO (VRE1261). The versions 
described in this paper are versions JAAJO 1 000000, 
JAAKO 1 000000, JAALO 1 000000 and JAAMO 1 000000. 

Results 

PFGE patterns and antimicrobial susceptibility of VREfm 
and VSEfm 

In an attempt to identify a putative ancestor of the VREfm 
outbreak strain, we examined consecutive blood isolates (n = 191, 
Collection B) of E. faecium from Karolinska University Hospital, 
Huddinge, during a four-year period (2006-2009) covering the 
time of appearance of the presumably first clinical isolate of the 



outbreak VREfm strain. A dominant PFGE pattern (n = 37; 26 
VSEfm and 1 1 VREfm) was indistinguishable or closely related to 
that of the VREfm outbreak strain (named SE-EfmB-0701). 
Representative isolates with this PFGE-pattern are shown in 
Figure 1. Among the blood isolates analysed in retrospect, this 
PFGE pattern was observed for the first time in a vancomycin 
susceptible isolate (VSEfm) from February 2007 and soon after in 
two more VSEfm. Three additional VREfm with identical PFGE 
pattern were detected during the autumn of 2007. The first known 
clinical isolate (VRE0651) was found in an abdominal infection in 
August 2007. In 2008 another eight VREfm and 14 VSEfm 
isolates with the same PFGE pattern were detected. During 
January until August 31'^ 2009 no more VREfm but still 10 
VSEfm of the SE-EfmB-0701 PFGE type were detected. 

Antimicrobial susceptibility testing of the 37 blood culture 
isolates revealed high-level gentamicin resistance (n= 14; 38%) as 
well as ampicillin (100%) and ciprofloxacin (100%) resistance and 
teicoplanin susceptibility (100%). The vanB gene was detected in 
all VREfm isolates. 

In collection A (n = 21), representing PFGE patterns as diverse 
as possible at the time of selection (autumn 2008), a total of 9 
subgroups (a to i) of the pattern SE-EfmB-0701 were found 
(Table 1). All PFGE patterns displayed a similarity >90%, thereby 
fulfilling the suggested definition of relatedness (>81%) according 
to Morrison et al. [38]. The isolates did however not group 
together consistently in relation to their geographical origin 
(Figure 1). 

The nineteen VREfm isolates in collection A had vancomycin 
MICs ranging from 8 to >256 mg/L and were susceptible to 
teicoplanin, consistent with the vanB2 genotype. All the isolates 
were resistant to ampicillin and ciprofloxacin, and seven displayed 
high-level resistance to gentamicin. The two VSE-isolates 
(VSE1036 and VSE1027) in collection A were susceptible to both 
vancomycin and teicoplanin, but expressed resistance to ampicillin 
and ciprofloxacin (Table 1). 

MLST and WGS analyses of VREfm and VSEfm 

Collection A (n = 21) were studied in greater detail by MLST 
and PGR for virulence genes (Table 1). Three isolates of the 
outbreak strain and one unrelated VREfm isolate were examined 
by WGS; VSE1036 (ST 192), VRE1044 (ST 192), VRE1261 
(ST192) and VRE0576 (ST17). In coUection A, all isolates 
including the two VSE belonged to ST 192 except one single 
locus (VRE0673, ST78) and two double locus (VRE0881 and 
VRE0576, ST17) variants of ST192. All isolates except VRE0881 
and VRE0576 belonged to the ST78 lineage. All isolates except 
VRE0576, a pre-outbreak vanB-positive ST 17 isolate from 2006 
had a related PFGE pattern as shown by Smal PFGE (Figure 1). 
The MLST- and PFGE results concurred moderately in showing 
relation, since the isolates with the most divergent PFGE patterns 
(VRE0673, VRE0776 and VRE0881) in two of three cases had 
deviating MLST profiles. 

The WGSs from the outbreak isolates VRE1044, VRE1261, 
and VSE 1036 as well as the unrelated vanB-positive pre-outbreak 
isolate VRE0576 were aligned against the chromosome of an E. 
faecium isolate from Australia which contains a vanB2 transposon 
(Aus0004) (Figure SI, vanB2 transposons are indicated by red 
triangles). The WGS of the three ST 192 outbreak isolates were 
naturally more homologous and shared several regions (highlight- 
ed by black triangles) that were not present in the genome 
sequences of either VRE0576 or Aus0004 (both ST 17) although it 
should be noted that the plasmid sequences of Aus0004 were not 
included in this comparison. The VRE0576 and Aus0004 
genomes showed many unique regions (white areas). Furthermore, 
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Lane Isolate no. County 



PFGE-type 




VRE0726 
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SE-EfmB-0701d 
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Figure 1. Dendrogram of Sma\ PFGE of the 21 isolates in collection A. The dendogram shows that 20 of the isolates are clonal (lane 1-20) 
and one is divergent (lane 21). The symbol ♦ in the dendrogram indicates a similarity of 90.6%. The PFGE-type nomenclature is based on the 
following: SE stands for Sweden, EfmB stands for E faecium with vanB, the number 07 represents year 2007 (the year the index was identified) and the 
last number is a serial number. The letter at the end describes which PFGE-subtype the isolate belongs to. 
doi:1 0.1 371 /journal.pone.01 03274.g001 
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VRE1044 and VRE1261 showed some unique regions (green 
triangles) that were not present in VSE1036. Further analyses of 
these regions suggested they belong to mobile genetic ele- 
ments. VRE 1044 and VRE 1261 showed only minor difFerences. 

Antimicrobial resistance determinants and virulence 
genes of VREfm and VSEfm 

All VREfm isolates in collection A harboured the vanB2 gene 
as an integral part of Tnl349/Tn3382 demonstrated by the 
mnX^-ORFC-PCR (Table 1). This link was also confirmed by 
WGSs in three VREfm isolates (Table S2). WGS data of VRE576 
revealed the same genetic organisation of the vanB2 transposon as 
well as 99% nucleotide (nt) identity to Tnl549. Interestingly, 
WGS data from VRE 1044 and VRE 1261 showed the same 
vanB2 transposon organisation and 99% nucleotide identity to 
Tn75^9, but also an additional 2588 bp inserted between nt 5014 
and 5015 of Tn75^9. We performed ICE^Wn Q8 PGR (Table 
SI) covering the putative insertion region in the remaining VRE 
isolates to disclose a potential unique insertion signature of the 
vanB transposon in SE-EfmB-0701 isolates. Presence of an 
approximately 2.6 kb insertion was confirmed in the vanB 
transposon of all the VREfm SE-EfmB-0701 isolates in collection 
A. The 2588 bp insert sequence is 89% identical to the region in 
Clostridium saccharolyticum-like KIO (GenBank Acc. 
No. FP929037) encoding a retron-type reverse transcriptase. In 
line with this, the 2588 bp sequence encodes a putative protein of 
610 amino acids (aa) with 99% identity to a putative reverse 
transcriptase /maturase from Faecalihacterium prausnitzii A2-165 



(GenBank Acc. No. EEU96266) and a putative group II intron- 
encoded protein LtrA (reverse transcriptase and RNA maturase) 
from Flavonif factor plautii ATCC 29863 (GenBank Acc. 
No. EHM54980). The putative protein further shows 43% 
identity to the group II intron 599 aa multifunctional protein 
LtrA in Lactococcus lactis (GenBank Acc. No. U50902). LtrA is 
known to have reverse transcriptase, RNA maturase and site- 
specific DNA endonuclease activity mediating intron splicing and 
mobility [39]. 

The gyrA and parC genes extracted from the WGSs revealed 
SNPs associated with ciprofloxacin resistance. Two mutation 
events in each gene were found (Table 1), and both aa 
combinations (GyrA Arg83, ParC Ile80 or GyrA Ile83, ParC 
Arg80) have been described previously in E. faecium isolates with 
ciprofloxacin MICs >16 mg/L [25,40,41]. Moreover, the tetra- 
cycline resistance determinant tetM and the macrolide resistance 
determinant ermB were also found in three and four WGS isolates, 
respectively (Table S2). 

All isolates of the PFGE type SE-EfmB-0701 harboured esp, 
sgrA, acm, scm,pilB, efaAfm, or/2010 and orf 2514. Moreover, 17 
of 20 isolates contained hyl. PGR data showed that the genes pilA, 
echA and orf903 occurred in six, one and one of 20 isolates, 
respectively. The ecbA and orf 903 genes were found in isolates 
with unique PFGE subtypes in this collection (VRE0776 and 
VRE0673) (Table 1). WGS data revealed that pilA (VRE1044 and 
VRE1261) and ecbA (VSE1036, VRE1044 and VRE1261) were 
present with a nucleotide match of 1672/1976 (85%) for pilA and 
2766/3173 (87%) for ecbA compared to reference sequences in 
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isolates El 162 and TX16, respectively. In addition, VSE1036, 
VRE1044, and VRE1261 carried a truncated version of pilA 
(219 bp) 100% identical to the reference sequences. Blastn search 
of the pilA and ecbA sequences within the NCBI shotgun sequence 
database revealed that these two genes had a notable SNP 
variation between different strains (0 20%). The original j??iL4 and 
ecbA primers used yielded no good matches in the WGS sequences 
which explain why we got few positive PGR products with these 
primers. All our collection A isolates were positive for these genes 
using new primers (pilA 2 and ecbA 2 Table SI) targeting 
conserved regions in the WGS and reference ecbA and jf?zZ/4 genes. 

Plasmid localisation of Tn7549-type transposon 

By BLAST alignment of the contigs from VRE1044 and 
VRE1261 containing the vanB2 transposon against the WGS data 
of the VSE1036 isolate, the exact AT-rich location of the vanB2 
transposon insertion site could be identified in contig00062 of 
VSE1036 (Figure 2 and blue triangle in Figure SI). The 
transposon insertion site was identical for VRE1044 and 
VRE1261 and corresponded with 100% identity to sequence in 
contig00062 of VSE1036 (Figure 2). 

To get a better picture of the plasmid content we performed 
Southern blot hybridisation of SI nuclease treated total DNAs 
separated by PFGE using specific replication (rep) gene probes 
(Table 1). WGS data gave additional information on plasmid 
sequences as well as resistance genes and Tnl349/Tn3382 (Table 
S2). Replication genes of r^jf?-classes ^^j^2/pRE255 ^^i^ii/pB825 ^^p\i/ 
pRUM and ^^/?pLGi were present in all four WGS, with some 
nucleotide differences in the pre-outbreak isolate (VRE0576) 
compared to the other three isolates. VRE0576 also differed from 
the other isolates by lacking a r^jf?i4/pRii -class gene present with 
identical nucleotide identity scores in the other WGSs, and by 
containing a putative ^^/>umque/pCiZ2"class gene absent in the other 
WGSs. Notably, VRE1044 and VRE1261 contained two rep2/ 
pRE25-sequences with SNP differences on different contigs, 
consistent with the observed two r^jf?2/pRE25 plasmids in the 
hybridisation analyses (Table 1 and S2). The VSE1036 isolate had 
a similar hybridisation pattern as VRE1044 and VRE1261, but 
did only have one 50-kb plasmid and one contig with ^^/>2/pRE25- 

Hybridisation results showed that all VREfm isolates harboured 
the vanB resistance gene on an approximately 70-kb ^^^17/pRUM 
replicon (shown for representative isolates VRE0726, VRE0734 
and VRE0881 in Figure 3 lanes 5, 8 and 11, and VRE0690, 
VRE0653, VRE0776 in Figure 4 lanes 5, 7 and 9). Further 
hybridisation of selected isolates revealed a dominant common 
plasmid content pattern for VRE0726, VRE0651, VRE0678, 
VRE1044, VRE1261, and VRE0650, and some minor differences 
in the other 1 1 characterised isolates (Table 1). 

Interestingly, the two VSE isolates clonally related to the 
VREfm outbreak strain (VSE1027 and VSE1036) differed from 
the common pattern by having a ^^j^iy/pRUM replicon 30 kb 



smaller than that of the VREfm isolates (Table 1), corresponding 
to the size of the vanB2-Tnl549/Tn5382 transposon. Notably, 
the pRUM rep A was found with 100% homology to the reference 
sequence in VRE0576 (GenBank Accession number 
JAAKOOOOOOOO), and present with 97% homology in VSE 1036 
(JAAJOOOOOOOO), VRE1044 JAALOOOOOOOO) and VRE1261 
(JAAMOOOOOOOO). A subsequent comparison of the r^jf?i7/pRUM 
replicon sequence in the three latter strains showed 100% 
homology. Thus, these plasmids share identical signature sequenc- 
es for pRUM rep A as well as for the Tnl549/Tn5382 transposon 
insertion site. 

Tn75^9/Tn535*2-related genes, axe-txe and the vanB2 gene 
were found within the same contig in the VRE 1 044 and VRE 1261 
isolates, thus supporting the SI nuclease PFGE hybridisation and 
PGR data showing linkage of vanB2-Tnl549/Tn5382 and axe- 
txe genes on the same replicon (representative isolates in Figure 3 
lanes 5, 8 and 11). WGS analyses of VSE1036 showed that the 
sequence corresponding to the transposon insertion site of 
VRE 1044 and VRE 1261 (Figure 2) mapped to the same 
contig00062 as axe-txe (Table S2) typically found on ?^^j^i7/pRUM 
replicons. Interestingly, the reference pRUM plasmid contains a 
putative relaxase gene and a mobilisation gene both consistently 
absent from all four WGSs. 

Transferability studies of the VanB determinant 

All but 2 (VRE0651 and VRE0673) out of 10 tested VREfm 
isolates produced transconjugants with a low frequency of 1 0 ^ to 
10 per donor, close to the detection limit, using E . faecium 64/ 
3 as a recipient (Table 2). Transconjugants were confirmed by 
PFGE using 5m^zI-digestion (Figure S2 and S3). E. faecium 
recipient BM4105-RF did not support conjugation within the 
detection limit. 

SI nuclease PFGE and subsequent hybridisation of the 64/3 
transconjugants revealed that all the transconjugants contained a 
vanB -rep ij/-p^UM plasmid of variable size (110-170 kb), larger 
than the original 70-kb donor plasmids (Figures 3, 4, S5 and S6). 
With the exception of VRE0690x64/3 (170 kb) (Figure 4 lane 4 
and Figure S5 lane 7) and VRE0776x64/3 (110 kb) (Figure 4 
lane 8 and Figure S5 lane 9), the transconjugants had a vanB- 
^^i^i7/pRUM-^^i^2/pRE25 Containing plasmid of about 140 kb 
(Figure 3 lanes 6, 9 and 12, Figure 4 lane 6, Figure S5 lanes 3, 
5 and 11 and Figure S6 lanes 5, 8 and 11) which corresponds to 
the joint size of the 70-kb r^j^2/pRE25 and the 70-kb vanB-repu/ 
pRUM of their donors. Conversely, the two donors VRE0690 and 
VRE0776 that gave transconjugants with different sized plasmids 
both lacked a copy of the 70-kb ^^/?2/pRE25 plasmid (Table 1). 
VRE0776 contain a 40-kb Tep2/pKE25 plasmid which joint with the 
70-kb vanB -rep 17 /p^uM could result in the 110-kb vanB-repu/ 
pRUM-^^j^2/pRE25 Seen in Figure S5 lane 9. In VRE0690 the 
170 kb plasmid of the transconjugant did not hybridise to rep2/ 
pRE25 suggesting that the 1 70-kb plasmid result from joining of the 



Left end Jnl549/5382 Right end 

VRE0576 ...atcacttgattttgttttacaattAAAATTTTAGGTTATCGCT6..y-33.8kb/.XTA2UU^TCCa^TAl^ 

II II I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

VSE1036 ...ctaaaatcagctaattttt ctttagaataaaaaaacatt... 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

VRE1261 ...ctaaaatcagctaatttttttctaAAAATTTTAGGTTATCGCTG..y-36.3kb/.XTAAAATCCCATATJU^TTTTctttagaataaaaaaa 

I I I I I 1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

VRE1044 ...ctaaaatcagctaatttttttctaAAAATTTTAGGTTATCGCTG.../~36.31cb/...CTAAAATCCCATATAATTTTctttagaataaaa 

Figure 2. Sequence comparison of the insertion regions of Tn 1549/5382, The figure shows the transposon insertion regions of VRE0576 
versus VRE1044 and VRE1261 and the corresponding region in VSE1036 (contig00062). 1^1549/1^5382 left and right end imperfect inverted repeats 
are shown in bold capital letters. Vertical lines indicate identical nucleotides. 
doi:1 0.1 371 /journal.pone.01 03274.g002 
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kb 
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145.5- 

97.0- 
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vanB 




^^Pl7/pRUM 

1 2 3 4 5 6 7 M 8 9 10 11 12 13 




axe-txe 



Figure 3. S1-nuclease PFGE and corresponding Southern hybridisations with repiy/pRUM/ i^^a?^ and axe-txe probes. These results 
illustrated transfer from donors VRE0726, VRE0734 and VRE0881 (lanes 5, 8 and 1 1) of a similar sized plasmid (approximately 140 kb) to 64/3 (lane 1) 
(1^* generation transconjugants shown in lanes 6, 9 and 12) which was subsequently retransferred to BM4105Str (lane 2) (2"^ generation 
transconjugants shown in lanes 7, 10 and 13) when using the 1^* generation transconjugants as donors. Lane 3 vanB positive control V583, lane 4 
''epi7/pRUM/ oxe-txe and vanB positive control £ faecium U37, lanes M low-range PFGE marker. 
doi:1 0.1 371 /journal. pone.01 03274.g003 



70-kb vanB-repij /p^uM with the 100-kb unknown replicon 
(Figure S5 lane 7). 

The transconjugants were tested for susceptibility to streptomy- 
cin in order to identify eligible donors for re transfer experiments. 
The first generation transconjugants originating from donors 



VRE0726 (ST 192), VRE0734 (ST 192) and VRE0881 (ST 17) 
were susceptible. The streptomycin resistance in the other 
transconjugants probably originate from co-transfer of streptomy- 
cin resistance with the plasmids of their clinical isolate donors 



M1 23456 789 M1 23456789 M1 23456789 




S1 nuclease PFGE 



vanB 



''^Pl7/pRUM 



Figure 4. S1-nuclease PFGE and corresponding Southern hybridisations with vanB and A'epi7/pRUM probes. Donors VRE0690, VRE0653 
and VRE0776 (lanes 5, 7 and 9) and their respective transconjugants (lanes 4, 6 and 8) illustrate transfer of different sized plasmids co-hybridising to 
vanB and repu/pRUM (circa 110-150 kb) into 64/3 (lane 1). Lane 2 vanB positive control V583, lane 3 repij/pR\jM and vanB positive control E. faecium 
U37, lane M low-range PFGE marker. 
doi:1 0.1 371 /journal.pone.01 03274.g004 
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Table 2. Transfer frequencies between donors and recipients after filter mating. 





Primary matings 




Transfer frequency 






Donor (PFGE 
subtype) (ST) 


Recipient 


Transconjugants/donor 


Transconjugants/recipient 




VRE0653 (a) (192) 


64/3 


6x10^ 


4x10^° 




VRE0651 (b) (192) 


64/3 


<6x10^^ 


<2x10^^ 




VRE0683 (c) (192) 


64/3 


8x10^ 


5x10 " 




VRE0734 (c) (192) 


64/3 


2x10 " 


2x10^^ 




VRE0726 (d) (192) 


64/3 


1x10^ 


3x10 " 




VRE0690 (e) (192) 


64/3 


2x10^ 


1x10^° 




VRE0688 (f) (192) 


64/3 


7x10^ 


6x10^ 




VRE0673 (g) (78) 


64/3 


<2x10^^ 


<1 xlO^^ 




VRE0776 (h) (192) 


64/3 


1x10^° 


5x10^ 




VRE0881 (i) (17) 


64/3 


2x10^ 


7x10 " 




Secondary matings 


VRE0726X64/3 




BM4105-Str 


7x10^ 


2x10^ 


VRE0734x64/3 




BM4105-Str 


1 xlO^ 


2x10^ 


VRE0881 X64/3 




BM4105-Str 


3x10"^ 


1 xlO^ 



^Isolates VRE0651 and VRE0673 did not support conjugation within given detection limits. 
doi:1 0.1 371/journal.pone.01 03274.t002 



(VRE0653, VRE0683, VRE0690, VRE0688 and VRE0776) that 
all showed high-level resistance to streptomycin. 

Retransfer and following SI nuclease PFGE and Southern 
hybridisations demonstrated that the vanB-rep u/^^uM-'^ep 2 /pKE25 
140-kb hybrid plasmids (Figure 3 and Figure S6) were stable in 
size and readily transferable with transfer rates of 10 "^-10 ^ 
transconjugants per donor (Table 2). Transconjugants were 
confirmed by PFGE using ^Wttl-digestion (Figure S4). 

Discussion 

A total of 872 VREfm cases were notified during the outbreak 
in 2007 to 2011, predominantly as faecal colonisation in elderly 
hospitalised patients with underlying diseases. Less than 10% of 
the VREfm isolates were recorded from blood, urine or wound 
samples [15]. 

Molecular characterisation of a representative subset of related 
outbreak isolates showed a strong predominance of ST 192, which 
is a single-locus variant of ST 7 8 and considered to be a high-risk 
genotype [42]. The MLST- and PFGE results were moderately 
congruent, as the isolates with less similar PFGE patterns and 
virulence profiles (VRE0673, VRE0776 and VRE0881) in two of 
three cases had a divergent MLST profile. Previous BAPS data 
have concluded that the ST78 and ST 17 lineages are located in 
BAPS group 2-1 and 3-3, respectively, and not closely related. 
This contrasts our PFGE results which grouped VRE0881 (ST 17 
lineage) in the SE-EfmB-0701 PFGE type (ST78 lineage). This 
observation as well as the carriage of a unique virulence gene by 
VRE0673 {orf903) and VRE0776 (only ecbA gene detected by the 
original PGR primers), could be explained by the lean inclusion 
criteria of relatedness by PFGE, thereby accepting not genetically 
related pulsotypes as related. The long collection time (around a 
year) and evidence of increased DNA banding pattern polymor- 
phism by E. faecium compared to other bacteria [38] was used as 
basis for our choice of inclusion criteria. However, the unique 
vanB transposon signature found in all the SE-EfmB-0701 isolates 
links VRE0881, VRE0776 and VRE0673 to the outbreak. 



The examined isolates expressed multidrug-resistance and 
harboured several specific genes associated with increased 
virulence. The vancomycin susceptible isolates (VSE1027 and 
VSE1036) from the start of the outbreak period were clonally 
related and exerted the same co-resistance- and virulence profile as 
the VREfm isolates (Figure 1 and Table 1). Moreover, plasmid 
profiling and WGS data (Table 1 and S2, Figure SI) also indicated 
close relatedness between SE-EfmB-0701 PFGE type isolates with 
minor differences in plasmid profile as well as clear differences in 
gene content compared to the pre-outbreak isolate. 

Our results support the notion that internalisation of the vanB 
transposon into the ^^/>i7/pRUM plasmid coincide with the 
successful spread of this high-risk strain. The ^^/>i7/pRUM replicon 
has previously been shown to harbour a segregation stability 
module encoded by a toxin-antitoxin cassette {axe-txe) which have 
been shown to support maintenance of linked antimicrobial 
resistance genes [43]. Notably, :r^j^i7/pRUM replicons with the axe- 
txe cassette have been shown to be present in a majority of CGI 7- 
like strains [22]. WGS analyses revealed that the SE-EfmB-0701 
pRUM replicons contained an axe-txe module with a 100% 
identity to the original pRUM axe-txe sequence. 

A repi7/p^uM-yG,nB2-Tnl349/Tn3382-axe-txe-p\^smid of ap- 
proximately 120-130 kb has previously been described in a 
polyclonal cluster of E. faecium from 2002-2004 in the Swedish 
county Orebro. This cluster originated from BAPS-group 3-3 
(ST 17, ST 18 and single locus variants of these) [18]. The Orebro 
isolates displayed a different PFGE pattern from SE-EfmB-0701 
[15] and had a vanB2 transposon without the unique signature 
found in SE-EfmB-0701. 

The repeatedly observed fusion between the t;(2n52-containg 
^^J^i7/pRUM plasmids and ^^/>2/pRE25 replicons during conjugation 
experiments is an interesting feature. The WGS data did not 
support the presence of the putative relaxase and mobilisation 
protein associated with the reference pRUM plasmid. This could 
explain the need for the vanB 2 -conta.ining repu/p^uM plasmids to 
fuse with a conjugation system from other intracellular sources in 
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order to be mobilised. Several studies have described mosaicism 
and/or recombination events between enterococcal plasmids 
[22,44] which may support enhanced host range or other 
functional benefits associated with several replicons in one 
plasmid. 

There is a noteworthy reservoir ofvanB in intestinal anaerobes, 
and introduction of vanB2-Tnl349/Tn3382 in enterococci from 
other co-habitants (mainly Gram-positive anaerobes) in the 
intestinal environment has been experimentally observed [45]. 
Howden et al. [46] tested the ecological impact by phylogenetic 
analysis of the transposons and their insertion sites, and showed 
that a diversification was likely due to a higher grade de novo 
VRE generation compared to cross-transmission between entero- 
coccal strains than previously believed. They also observed an 
increasing incidence of nosocomial VRE infections despite 
engagement of control interventions to limit transmission between 
patients. Based on the extensive PFGE and selective MLST 
analyses in this study, it is highly probable that the closely related 
VSEfm ST 192 strain was a successful hospital coloniser in Sweden 
already in 2007. Acquisition of the vanB2 transposon by the 
VSEfm ST 192 outbreak strain is the most likely hypothesis on how 
vancomycin resistance appeared in this strain. The WGS data 
strongly support that a vanB2 transposon with unique signature 
was inserted within a r^jf?i7/pRUM-plasmid with an unique pRUM 
repA sequence signature in a strain already present, causing a 
parallel evolution between VSEfm clones without, and VREfm 
clones with the rep 17 /p^uM-'^^(^nB2-Tnl 349 /Tn3382-axe-txe ar- 
rangement. 

Suggested clearance time for VRE faecal colonisation is 
estimated to be 4 years [47], but others have suggested that 
environmentally adapted VRE are capable of inhabiting the 
intestines in small numbers for even longer [48,49]. 

The selective enriched broth used in some laboratories in 
Sweden before January 2009 with a vancomycin concentration of 
32 mg/L was not suitable for vanB-type resistance screening. To 
address this problem, microbiological laboratories were then 
advised to reduce the concentration to 4 mg/L [15]. However, 
vanB-type VRE may have even lower MICs [17,50]. Importantly, 
the EUCAST disk diffusion test used by most laboratories in 
Sweden as the phenotypic vancomycin-susceptibility test method 
rely on observing the zone edge quality for identification of low 
level vanB-typG resistance. This introduces observer experience as 
a variable [51]. 

In conclusion, the molecular analyses revealed that the E. 
faecium outbreak strain belonged to the high-risk genetic lineage 
of ST 192. The strain was resistant to several commonly used 
antibiotics and harboured several virulence genes. A successful 
^^j^i7/pRUM~pl3^smid containing a vanB transposon with a unique 
genetic signature originating from other intestinal bacterial species 
was present in all the VREfm isolates related to the outbreak 
strain. The ^^j^iy/pRUM plasmid harboured a toxin- antitoxin 
module supporting plasmid maintenance. In addition the repu/ 
pRUM replicon can easily join with conjugative genetic elements 
supporting spread to other high-risk E. faecium clones. The 
current phenotypic screening methods might hamper efforts in 
limiting VREfm spread, as low-MIC vanB -type VRE might go 
undetected [51]. 

Supporting Information 

Figure SI Alignment of genome sequences from VRE 1044 (row 
2), VRE 1261 (row 3), VSE1036 (row 4) and VRE0576 (row 5) 
against the reference genome Aus0004 NG_0 17022 (row 1). The 
similarity plot indicates average similarity for each region. 



Coloured blocks indicate regions of sequence homology in the 
genomes and white areas indicate regions with low sequence 
homology. Red triangles indicate contigs containing Tnl349/ 
Tn3382 (VRE 1044 contigs 00036 and 00041, VRE 1261 contigs 
00049 and 00044/VRE0576 contig 00004). The blue triangle 
indicates the transposon Tnl349/Tn3382 insertion region in 
contig 00062 of VSE1036. This contig also contains the axe-txe 
genes typically found on ^^/>i7/pRUM replicons. Black (VRE 1044, 
VRE1261 and VSE1036) and green (VRE1044 and VRE1261) 
triangles highlight regions with contigs or partial contigs found in 
the outbreak isolates but not in VRE576 or Aus0004. 
(PDF) 

Figure S2 Smal PFGE of first generation transconjugants (TG) 
(lanes 5, 7, 9, 11) showing divergent band patterns compared with 
the clinical isolate donors (lanes 6, 8, 10 and 12) and similar 
pattern with recipient 64/3 (lane 2). Lanes 1 and 13 low-range 
PFGE marker, lane 3 vanB positive control E. faecalis V583, lane 
4 ^^j^i7/pRUM positive control E. faecium U37, lanes 5 and 6 TC 
and donor VRE0726, lanes 7 and 8 TC and donor VRE0734, 
lanes 9 and 10 TC and donor VRE0683, lanes 1 1 and 12 TC and 
donor VRE0688. 
(PDF) 

Figure S3 Smal PFGE of first generation transconjugants (TC) 
(lanes 5, 7, 9, 11) showing divergent band patterns compared with 
the clinical isolate donors (lanes 6, 8, 10 and 12) and similar 
pattern with recipient 64/3 (lane 2). Lanes 1 and 13 low-range 
PFGE marker, lane 3 vanB positive control E. faecalis V583, lane 
4 ^^j^i7/pRUM positive control E. faecium U37, lanes 5 and 6 TC 
and donor VRE0690, lanes 7 and 8 TC and donor VRE0653, 
lanes 9 and 10 TC and donor VRE0776, lanes 1 1 and 12 TC and 
donor VRE0881. 
(PDF) 

Figure S4 Smal PFGE of second generation transconjugants 
(TCs) (lanes 6-8, 11-13 and 15—17) showing divergent band 
patterns compared with the first generation transconjugant donors 
(lanes 5, 10 and 14) and similar pattern with recipient BM4105-Str 
(lane 4). Lanes 1, 9 and 18 low-range PFGE marker, lane 2 vanB 
positive control E. faecalis V583, lane 3 ^^/>i7/pRUM positive 
control E. faecium U37, lane 5 donor VRE0726 x64/3, lanes 6-8 
TCs VRE0726x64/3xBM4105-Str, lane 10 donor 
VRE0734X64/3, lanes 11-13 TCs VRE0734x64/3xBM4105- 
Str, lane 14 donor VRE0881 x64/3, lanes 15-17 TCs 
VRE0881 x64/3xBM4105-Str. 
(PDF) 

Figure S5 SI -nuclease PFGE and corresponding Southern 

hybridisations with Tep2/pKE25 ^^/>i7/pRUM probes showing 
co-hybridisation in first generation transconjugants (lanes 3, 5, 9 
and 11). Lanes 1 and 12 low-range PFGE marker, lanes 2 and 3 
donor and TC VRE0683, lanes 4 and 5 donor and TC VRE0688, 
lanes 6 and 7 donor and TC VRE0690, lanes 8 and 9 donor and 
TC VRE0776, lanes 10 and 11 donor and TC VRE0653. 
(PDF) 

Figure S6 SI -nuclease PFGE and corresponding Southern 

hybridisations with Tep2/pKE25 ^^/>i7/pRUM probes showing 
co-hybridisation in first (lane 5, 8 and 1 1) and second generation 
transconjugants (lane 6, 9 and 12). Lane M low-range PFGE 
marker, lane 1 Tep u/^^um and Tep2/pKE25 positive control E. 
faecium U37, lane 2 recipient 64/3, lane 3 VRE1044, lane 4 
VRE0726, lane 5 VRE0726 x64/3, lane 6 VRE0726x64/ 
3xBM4105-Str, lane 7 VRE0734, lane 8 VRE0734x64/3, lane 
9 VRE0734x64/3xBM4105-Str, lane 10 VRE0881, lane 11 
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VRE0881 X64/3, lane 12 VRE0881 x64/3xBM4105-Str, lane 13 

recipient BM4105-Str. 

(PDF) 

Table SI Primers used in this article. 
(PDF) 

Table S2 Plasmid replication, resistance, toxin-antitoxin system 
and conjugative transposon genes found in the WGSs of the pre- 
outbreak isolate VRE576 and the three outbreak isolates 
VSE1036, VRE1044 and VRE126L Gene identity refers to the 
reference sequence. 
(PDF) 
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